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On 

Abstract. We study entanglement generation via particle transport across a one- 
dimensional system described by the Bose-Hubbard Hamiltonian. We analyze how the 
^ , competition between interactions and tunneling affects transport properties and the 

creation of entanglement in the occupation number basis. Alternatively, we propose 
to use spatially delocalized quantum bits, where a quantum bit is defined by the 
\ presence of a particle either in a site or in the adjacent one. Our results can serve as 

a guidance for future experiments to characterize entanglement of ultracold gases in 
. one-dimensional optical lattices. 

' 1- Introduction 

a 

q-i The generation of entanglement between distant nodes of a quantum network has 

profound implications for quantum computation and information pQ, and has triggered 
a remarkable effort into the study of entanglement generation and quantum state 
transport between distant lattice sites (see, e.g., Refs. [21 EH BJ E] for discussions of 
transport in spin chains; entanglement dynamics in such systems is, e.g., studied in 
OtZHH]). Furthermore, the recent experimental achievements in loading either bosonic 
and/or fermionic ultracold atomic gases into optical lattices, which permits to reproduce 
very accurately several spin Hamiltonians, have spurred an enormous interest in lattice 
systems (see [9] and references therein). 

Here we focus on one of the simplest but yet non-trivial lattice model, the so-called 
one-dimensional (ID) Bose-Hubbard Model (BHM) [TUl [TT]. It describes a system of 
spinless bosons with (repulsive) on-site interaction, which can hop (tunnel) between 
adjacent sites of a ID lattice. In general, this genuinely many-body model cannot be 
reduced to an effective non-interacting one. As proposed in [12] and later experimentally 
demonstrated [13], it can be realized by confining ultracold atoms in an optical lattice. 
In this article we study transport properties in this system. Here, by transport we refer 
to the dynamics obtained when an extra particle is loaded onto a system that previously 
was cooled to its ground state. In some limiting cases, transport can be described in the 
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language of continuous time quantum walks (a quantum analogue of classical random 
walks) and closed analytical expressions can be found [H]. When competition between 
hoping and on-site interactions arises, the transport properties as well as the generation 
of entanglement between distant locations of the lattice is substantially modified. 

The article is structured as follows: in the present section we review the essential 
properties of the Bose-Hubbard model. As one limiting case we identify transport 
across the lattice with continuous time quantum walks. We discuss quantification and 
characterization of entanglement for a single particle propagating in a ID lattice. In 
Sec. |2j we analyze the generation of entanglement between distant lattice sites when an 
extra particle is loaded on top of the ground state. In Sec. El we discuss entanglement 
in the so-called Spatially Delocalized Qubit (SDQ) basis and investigate the role of 
interactions within this scheme. Finally, we summarize briefly in Sec. HJ 



1.1. The Bose-Hubbard Hamiltonian 

The Bose-Hubbard Hamiltonian for a ID lattice of M sites (with open-boundary 
conditions) has the form 

M-l M jj M 

i=l i=l ^ i=l 

where dj and a\ are the bosonic annihilation and creation operators for a particle on the 
z'-th lattice site, hi = djd; is the corresponding bosonic number operator, and e« accounts 
for the single-particle on-site energy. The Bose-Hubbard model assumes only nearest- 
neighbour tunneling with constant amplitude J and pairwise interaction between bosons 
on the same site leading to an energy shift U. 

A particularly clean realization of such a Hamiltonian is obtained by trapping 
neutral bosons via the dipole force in an optical lattice [12]. By taking the trapping 
sufficiently tight in two directions (say y and z), an effective one-dimensional system 
can be realized. The corresponding Hamiltonian in second quantization notation reads 

U OL = J dxtjj j (x) + Vosin 2 (7rx/c/)^ ip(x) +g J dxft '(x)tp ] ' (x)ijj(x)i/j(x) 

= j dxj)\x)H b Jj){x) +74*. (2) 

The optical lattice is completely characterized by its depth Vq, that can be controlled 
through the laser intensity, and by the wave number k = 2n/\ where d = A/2 is 
the lattice periodicity. The effective ID interaction strength g = 2irTia s uj t (where u t 
is the transversal frequency of the trap assumed equal in the y, z directions and a s is 
the atomic scattering length) can be changed either by modifying the confinement of 
the atoms in the two orthogonal directions or, alternatively, via a Feshbach resonance, 
which even allows to change the sign of a s jT5j. For periodic boundary conditions 
(or large enough lattices) the bosonic operators can be expanded in terms of Bloch 
functions. In the low temperature regime and with typical bosonic interaction strengths, 
excitations to higher bands can be neglected for sufficiently deep lattices. The dynamics 



Transport and Entanglement Generation in the Bose-Hubbard Model 



3 



is then restricted to the lowest Bloch band and the field operators can be expanded 
in terms of single-particle Wannier functions localized at each lattice site Xi = di as 
w(x — Xi) = (x\wi). The bosonic creation operators a\ are now defined via \w/) = a\\Q), 
\n) being the vacuum. Eq. ([1]) is derived from Eq. ([2]) by keeping only nearest-neighbor 
hopping and restricting the interactions between bosons to a contact (zero-range) 
potential. Under these approximations, the tunneling amplitude between adjacent 
sites reads J = (wi\Hf TCC \w i+ i) , and the on-site boson-boson interaction is given by 
U = g J dx\w(x)\ A . 

The Bose-Hubbard ground state with filling factor n is obtained by minimizing 
(Hbr—h Yji n i)i where the chemical potential /i fixes the total number of particles. In the 
limit U/J^O (strictly speaking, in the thermodynamic limit for U / J < (U/ J) c ~ 3.44 
for n = 1 [16]) it is energetically favorable to spread each particle over the whole lattice. 
For periodic boundary conditions, a ground state with filling factor n can be explicitly 
written as 

Mn 

|n>, (3) 

where Mn is the total number of particles. This superfluid (SF) state is characterized 
by large fluctuations of the on-site number of particles, divergent correlation length, 
and a vanishing gap. In the opposite limit, for U/J — > oo, the ground state is a Mott 
Insulator (MI) state, i.e., a product state with well-defined number n of atoms per site 

-i M 

iv£s,Mi> = ^n(<fi^ w 

where n — 1 < [i/U < n. The MI state has finite correlation length and a gapped 
spectrum. Increasing U/J from to oo for integer filling factor (at T = and in the 
infinite system), the lD-Bose-Hubbard model undergoes a quantum phase transition 
(which corresponds to a Kosterliz-Thouless phase transition). There is also a generic 
phase transition which is crossed when the value of U/J is fixed and the chemical 
potential \i changes. In this case, the number of particles is not conserved and the 
behavior of the ground state near the phase transition simply corresponds to a weakly 
interacting condensate with a non-integer filling factor. See (9] and references therein 
for a detailed review of the properties of the Bose-Hubbard model. 

1.2. Continuous time quantum walks 

For the simplest transport case in the Bose-Hubbard model, namely, for a single boson 
placed in an otherwise empty lattice, the dynamics is equivalent to the one of a free 
particle moving in (finite) discretized one-dimensional space. Recently, such a model 
and its generalizations to more complex underlying graphs have been studied intensively 
in the context of Continuous Time Quantum Walks (CTQWs). Quantum walks, either 
continuous or discrete, have been proposed as quantum versions of classical random 
walks and analyzed with the aim of constructing new types of algorithms. They have 
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also been studied, e.g., in relation to decoherence properties of lattice systems. See [T7] 
(and references therein) for an excellent review of the topic. 

The definition of a CTQW is closely related to classical continuous time random 
walks [TBI • Let us consider the classical situation of a "particle" which can move on 
a set of vertices. The probability to jump from vertex % to another vertex j per unit 
time is denoted as J^-, with > if both vertices are connected and = otherwise. 
To conserve the total probability we demand J ti = — Jij- If Pi(t) is the probability 
of being at time t at vertex i, then 

^ = E<W*)- (5) 

Given a quantum state in the Hilbert space H. — % = 1 . . . M} spanned by the 
vertices, the similarity between Eq. (EJ) and the Schrodinger equation for the amplitudes 

<■#(*)), 

suggests to define a quantum analogue of the classical random walk by identifying Jy 
with the Hamiltonian matrix elements: Hij = (i\H\j) = J^. For vertices arranged on a 
finite line with M sites and constant nearest-neighbour transition probabilities = J 
in the graph, equation ([H]) is equivalent to the Schrodinger equation of a single particle 
evolving under the Bose-Hubbard Hamiltonian ([TJ). To generate entanglement in an 
effective way, transport should be symmetric, thus we impose an odd number M of sites. 
Placing the particle in the middle of the chain, i.e., preparing \ip(t = 0)) = a^ M+1 y 2 \Q), 
dynamics after a time t leads to the state 

M 

WOHE^WfiJin), (7) 



where the coefficients Ci(t) are given by [21 [3] 



Ci (t) = ^^y 




sm — sm 



7rki 
M + l 



cxp 



2i Jt cos 



kn 



M + l 



(8) 



The probability distribution pi(t) = |cj(t)| 2 , see Fig. [TJ is symmetric with respect to 
the point (M + l)/2. For times t such that pi(t) = pnit) 1, its standard deviation 
A = \J^2iPii 2 — (JZiPii) 2 grows linearly in time: A cx t [T7] . This is in strong contrast 
to a classical ID random walk, where A oc \/t. 



1.3. Characterizing entanglement 

Here we first discuss shortly how to characterize entanglement distributed in a CTQW. 
We will later generalize this discussion to the transport of a defect, namely an extra 
particle, loaded on top of the ground state of the Bose-Hubbard Hamiltonian. Since 
particles in the CTQW as well as in the spinless Bose-Hubbard model have no internal 
degrees of freedom, the notion of quantum bit (or quantum dit) and its extensions to 
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Figure 1. (a) Probability distribution |ci(i)| 2 versus time for a CTQW starting at 
site (M + l)/2 in a chain of M = 41 sites; (b) probability distribution at time t = 5 J. 



entangled states and entanglement measures have to be redefined. S. Bose has discussed 
the distribution of entanglement in the CTQW for a single particle as well as for N non- 
interacting indistinguishable bosons [Hj (corresponding to N [symmetrized] versions 
of the quantum walk), quantifying entanglement between the two outer lattice sites 
(say 1, M) in the occupation number basis. To this aim, the reduced density matrix 
of sites 1 and M is expanded in the basis {In^nju)}, where rij corresponds to the 
number of particles at site i. Within this basis, usual entanglement measures can be 
applied. Quantifying entanglement in this "second quantized" formalism was introduced 
by P. Zanardi [19] and subsequently intensively discussed in the literature [20ll2~Tl 122] [23] . 
For just a single particle, mapping its presence (absence) on site % to spin-up (spin- 
down) reduces the propagation of a single boson in the lattice to the dynamics of a 
single, initially localized, excitation in the XY spin-chain [2j [3] - This analogy between 
"space" and "spin-1/2" entanglement is, however, limited: a spin at site i can be in any 
superposition of "up" and "down"; |up); + |down);, while |n, = 0) + |n, = 1) has no 
physical meaning for massive particles. Also, in the occupation number basis, operators 
must correspond to a direct sum of operators acting in sectors with fixed number of 
particles. For these reasons, it is not clear whether entanglement in the occupation 
number basis can be identified at all with usual "spin" -entanglement. In particular, 
there are no protocols that convert a (highly) entangled state in the occupation number 
basis to a "spin" -entangled state containing the same amount of entanglement!]]. 

Despite its practical drawbacks, we investigate first the generation of entanglement 
in the occupation number basis, thereby merely using it as a tool to characterize 
transport in the system. For a single particle (i.e., the CTQW), the reduced density 
matrix of sites 1 and M is always of the form 

Put® = P (0) pP> , p(°) = l-2 Pl (t), pM = 2 Pl (t)\* + )(* + \, (9) 

where the upper indexes correspond to the total particle number, and |^/ + ) = (|1,0) + 
|0, l))/v / 2. To later be able to generalize to states with more than one particle, we 

X In [TJ] a method to convert "space" to "spin" entanglement is given for the coupled chain. It 
requires interaction between spins and final measurements over intermediate sites. Thus it is not a 
local conversion scheme. Furthermore, it allows to extract only up to one ebit. In [20], a local protocol 
is presented to extract an ebit from (|1, 0) + 10, l))/v / 2, which however needs a sink/source for particles. 
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measure the entanglement through the logarithmic negativity LN [21]. From Eq. 
LN(p 1A/ (t)) =log 2 ||p[ M (t)|| 1 



log 2 



2pi(t) + 



I a + (1 - 2 Pl (t))^P , a - (1 - 2 Pl (t))^ 



(10) 



2 V 2 

where a = 1 - 4pi(i) + 6pi(t) 2 , = 1- A Pl (t) + 8pi(t) 2 . Since we are considering open 
boundary conditions, the logarithmic negativity presents several local maxima in time 
(as well as periodic revivals) due to multiple reflections. In order to characterize the 
CTQW and subsequently the Bose-Hubbard model through the generated entanglement, 
we consider only its first maximum. 

2. Transport properties of the Bose-Hubbard model: from the Mott 
Insulator to the superfluid phase 

2.1. Transport of an additional particle on a Bose-Hubbard ground state background 

We now generalize the previous discussion transport of a single particle in an empty 
lattice system to an extra particle propagating on a background given by a Bose-Hubbard 
ground state. Let us start by discussing the two extreme cases (i) U/J —> oo (Mott 
insulator phase) and (ii) U/J — > (superfluid phase). 

(i) Mott Insulator phase. Adding an extra-particle at site (M + l)/2 to the MI 
ground state with an integer filling factor n leads to the initial state 

m = o» = %afii« s ,M,> = %afi (*)" mm du 

For U I J large, during time evolution the system will remain in the subspace spanned by 
states {\i n ) = a{|V>GS,Mi)/V« + 1, i = 1...M}. Noticing that (-Jal +l a i )a\\'ijjQ SM1 ) = 
— J(n+ 1)oJ+i|'0gs mi)' we fi n( i that the effective Hamiltonian, up to corrections of order 
(U/J)~\ reads 

H eS = -J{n +1)E [\in){(i + l)n| + \{i + l)n){in\] • (12) 

i 

Thus, \ip(t)) = c r(^) I*")' with c™ = Ci((n + l)t), i.e., from bosonic enhancement the 
propagation is n + 1 times faster than in the pure CTQW, while the magnitude of the 
distributed entanglement is as before. 

(ii) Superfluid phase. In the opposite limiting case, for U/J = 0, approximating 
the ground state of the system with open boundaries by the one for periodic boundary 
conditions, the initial state reads 



t r -i M ~' Mn 



i m 



1-0(0)) = ^m+^IV^sf) = a^== 

V( Mn ) ! 

where a is a normalization constant. Evolution of this state leads to 

a M 



|n>, (is) 



l^(*)> = -7==E c *(*) a IlV'Ss,sP> J (14) 
'(Mn)\ i=i 
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lU ,j;,a]] = — J(a] +1 + a]_x). Thus, the extra particle added on 
top of the superfluid ground state propagates as a single particle in an otherwise empty 
lattice. 

The entanglement between the two outer sites however is strikingly different for the 
Mott insulator and the superfluid state, as we will discuss now. Let us first consider 
entanglement between the outer sites just for the ground state, without any additional 
particle. In the insulator case, we have 

1 M i i M-l 



I^GS.Ml) 



IK 



lt\n 



|n> = ^(* 

n! 



tftt \n 

m) 



(Vn!) 



n 



(15) 
(16) 



, ' n '-y- i=i 

= |0Ml)lM ® |^Ml)2...M- 

The reduced density matrix of sites 1, M then reads 

PMi.fi = fr2...A/-i(|^GS,Mi)(^GS,Mil) = \<ha) {<ha\ = \n,fi)(n,n\, (17) 

where the latter expression is in the occupation number basis, pum is a non-entangled 
pure state. On the other hand, for the superfluid ground state and M > 2, 



IV'gs.sf) 



(MS) 



M 



M 
a\ 



(18) 



M-l 



fc=o V 

Mn 

= ETfcl^F)^ ® |^SF)2...A/ 



V2 J y/(Mfi-k)\ \VM-2 p 2 




k=0 



with 



Ik 



fMn\ 2 k (M - 2) M "" fc 
\\k ) M Mn ' 



(19) 



(20) 



(21) 



The reduced density matrix reads 

Mn 

PSF,n = E7fc|0SF)(0S fc F l- 
k=0 

Except for the trivial cases Mn = or M = 2, the state Psf,h is mixed, has a direct 
sum structure in the occupation number basis, and is entangled. As Fig. [2] shows, 
entanglement in the ground state grows (weakly) with the total number of particles Mn 
in the system (keeping the number of sites M fixed). This happens despite the fact 
that the purity of the reduced state decreases at the same time, meaning that the outer 
sites also become entangled to the inner part of the chain. The reason is that adding 
more particles corresponds to adding more degrees of freedom which can be entangled 
in the occupation number basis. For this particular case, this leads to an increase in 
entanglement. 
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Figure 2. Entanglement (as measured by the logarithmic negativity LN) between the 
outer sites for the superfluid ground state (squares) and the superfluid ground state 
with an extra particle using Cj = (5a + <5,;m)/v / 2 (filled circles) as discussed in the text. 
The number of sites is M = 20. The purity tr[(/5sp) 2 ] of the reduced density matrix is 
plotted as a dashed line. For the insulator case (with integer n), always LN(p' Ml n ) = 
for the ground state and LN(p' MI n ) = 1 for the ground state with an additional particle 
and c, as above. 



In order to discuss the entanglement generated from an extra particle on top of the 
ground state, let us consider the simplified situation q = (Sn + 5jjvf)/\/2- For the Mott 
case, 

PMi,n = \ (I™ + 1, n) + \n, n + 1)) «n + 1, n| + (n, n + 1|) , (22) 

which is a pure state with logarithmic negativity LN(p' Mlfi ) = 1, independently of n. 
Note that the entanglement is completely contained in the sector of 2n + 1 particles 
shared between the outer sites. This is different for the superfluid phase. Here the extra 
particle on top of the ground state with filling factor n leads to 

Mn+l 

PsF, fi = E K)Vs fe F>(</>sU (23) 

k=0 

with 7q = 0, and 

K) 2 = ^^ for k>0. (24) 

Clearly, for an empty lattice, n — 0, the choice of q leads to LN(/3g Fn=0 ) = 1. For 
n > 0, the logarithmic negativity p' SF decreases (Fig. [2]). Again, as n grows the number 
of degrees of freedom which potentially can be entangled increases also. Still, in this 
case the two outer sites are less entangled due to the smaller purity of the reduced 
density matrix (the outer sites are more entangled to the inner part of the chain). It 
might seem counter-intuitive that 

LN(p' SFifi ) < LN(p' SFifi=0 ), (25) 

despite the additivity property of the logarithmic negativity [21]. Here we should again 
remark, that the occupation number basis which is used to quantify entanglement does 
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Figure 3. Entanglement between the two outer sites of the lattice, measured by 
the logarithmic negativity LN(/3im), versus U/J for a chain of M = 41 sites (after 
inserting an extra particle at site at (M + l)/2). The chemical potential /i is adapted 
to have a ground state with mean number of particles per site n = 1 (jx = U /(2J) for 
U 3> J). The data is obtained through MPS simulations truncating the basis to states 
of up to 5 particles per site and using D < 20. 



not reflect the tensor product structure of individual particles. In fact, the bosonic 
particles themselves live in a symmetrized subspace which does not have the structure 
of a full tensor product (this indeed is the reason why for bosonic particles new measures 
of entanglement have to be introduced). 

Let us finally note that if the systems dynamics is governed by the Bose-Hubbard 
Hamiltonian, then the logarithmic negativity will always be below the values of Fig. [2} 
as in this case the Cj(t)'s are given by Eq. (jSJ). 

To study how the competition between interactions and tunneling affects the 
generation of entanglement when adding an extra particle, we numerically calculate 
time evolution of the initial state for a wide range of parameters, 6 < U / J < 40. We use 
standard numerical MPS algorithms [25]. As we limit the number of particles per site to 
5 in our simulations, we cannot study states well inside the superfluid regime, where for 
n = 1 up to M particles per site have to be taken into account. Still, already for values 
of U/J < 15 drastic changes due to the possibility of tunneling are manifested in the 
entanglement between the outer lattice sites. To calculate \ip(t)) = exp(— ifzBH^I^GS 1 )' 
we first obtain the ground state of the Bose-Hubbard Hamiltonian for different values 
of the parameter U/J. We choose the chemical potential as \ij J = U/(2J) in the case 
of large U/J. For small values of the ratio U/J, we adjust fi appropriately, in order to 
obtain the ground state I^gs 1 ) w hh filling factor n = Jjj. Having the ground state, we 
add a particle to the system and obtain the dynamics through a time-dependent MPS 
simulations. From the MPS state, the reduced density matrix po/ of the outer sites 
can be extracted efficiently, such that we can compute the logarithmic negativity as a 
function of time. In Fig. [3] we display the value of LN(pi m) at its first maximum as 
a function of U/J. In Fig. [4] we analyze time evolution in detail for two cases: well in 
the MI phase (U/J = 40) and close to the SF phase (U/J = 6). Figs. H (a,c) show the 

§ In a one dimensional system the lobes of the Mott insulating phase are much stronger deformed 
than in two or three dimensions 1261 . 
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Figure 4. Left column: mean ocupation number rii = (hi) of site i as a function of 
time for a chain of M = 41 lattice sites ((a) U/J = 40; (c) U/J = 6). The system is 
initially cooled to its ground state with filling factor n = 1. At t — 0, a further particle 
is loaded at site i = 20. Right column: logarithmic negativity LN(p2i-fc,2i+fc) of the 
reduced density matrix of sites 21 — k and 21 + k, k = 1, . . . , 20 versus time for (b) 
U/J = 40, (d) U/J = 6. 

propagation of the excitation, i.e., the mean occupation number (n/) versus time. The 
propagation in the two cases is very similar (as it is visible from the figures, the evolution 
is slower for Uj J — 6, though not by a factor of 2 as it would be the case for Uj J — 0). 
The propagation of entanglement, visualized in Figs. H] (b,d) through the logarithmic 
negativities of the reduced density matrices of sites (M — l)/2±fc, is different in the two 
cases. Clearly, the efficient generation of entanglement in the occupation number basis 
requires a Mott insulator background. As we demonstrated, entanglement generation is 
much less efficient if tunneling becomes comparable to or larger than interactions. 

A strategy to increase entanglement in a system of ideal (interactionless) bosons 
in the occupation number basis consists in loading several bosons in a given site of an 
otherwise empty lattice [H] . For interacting particles, however, in the limit U ^> J, 
entanglement does not increase if all N particles are initially located at the same site. 
In this case, tunneling of a single particle alone is strongly suppressed and atoms tunnel 
together. Treating the tunneling term perturbatively, it can be seen that the evolution 
is slower by a factor J /(UN). 

3. Creating entangled "spatially delocalized quantum bits" 

In this section we discuss an alternative way of defining a quantum bit in a lattice filled 
with spinless bosons. This definition does not rely on the occupation number basis, 
leading to a notion of entanglement which is physically more sound. We use the concept 
of "spatially delocalized quantum bits" (SDQs), in which the binary alternative consists 
in having an atom either in one side or in the adjacent one. 

To be specific, we use a single particle shared between two (adjacent) sites i, i + 1 
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Figure 5. Probability p for the successful projection onto the subspace of one particle 
per spatially delocalized qubit (straight blue line) , logarithmic negativity LN (dashed 
green line) between the qubits after successful projection, and the product p LN 
(dashed red line) for a chain of 24 sites with (a) no interaction U/J = and (b) 
strong interaction U/ J = 20. Results from MPS simulations with D = 18. 



of the chain to define a quantum bit by identifying |0)sdq = \ n i — lj^j+i = 0) and 
|1)sdq = \ n i — 0,n i+1 = 1). Though in this definition again a qubit is defined via 
absence or presence of a particle on a lattice site, now one particle is necessary per qubit. 
In this way, for any unitary transformation the number of particles is conserved locally. 
For such "charge" or "spatially delocalized" qubits, implementations of quantum gates 
have been proposed for bosonic atoms in optical lattices |27j . as well as for electrons in 
quantum dots [28] or photons in photonic crystals |29j . 

To entangle two such qubits at the ends of a ID chain of even length M, we place 
two particles at sites M/2 and M/2 — 1 (in the middle of the chain) and let the system 
evolve. As before, we start by considering an otherwise empty chain, i.e., the quantum 
walk with two particles: \ip(t = 0)) = ^m/2^m/2+i\^) ■ For a spatially delocalized qubit 
defined through a particle shared between sites we introduce the corresponding 
projection operators = |a)sDQ(«|- For instance, 

Pi i,j) = |0) sdq (0| = k = hnj = 0)(n i = l, nj = 0|. (26) 
p^j) _|_ p(hj) ^ us projects onto the subspace having a single particle shared between 
sites i and j. If two qubits are defined on sites (1,2) and (M — 1,M) respectively, we 
obtain the density matrix for the two spatially delocalized qubits as 

pX'At) = p£' 2) ® p?-^m))m)\p2» ® p^~ 1M) - ( 27 ) 

The probability p of a successful projection onto the subspace of one particle per SDQ 
is given by p = trp sdq < 1, and the entanglement between the two spatially delocalized 
qubits in case of a successful projection is measured by the logarithmic negativity of 
the correctly normalized state LN(p sdq /p). The probability p, the logarithmic negativity 
LN, and the probabilistic entanglement p LN are plotted in Fig. [5] (a) for the case of no 
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Figure 6. Generation of entanglement between two spatially delocalized qubits 
formed from sites 1, 2 and 3, 4, respectively. Left column: U/J = 0, right column: 
U/J = 20. (a,d) show the mean site occupation (rii), (b,e) the populations 
SDQ(a/3|/S SD< ^(i)/p|a/3)sDQ, and (c,f) give the probability p of successful projection 
into the subspace of one particle per spatially delocalized qubit and the entanglement 
measured via the logarithmic negativity LN. 



interaction between the two bosons (U/J = 0) and for strong interaction between them 
(17/ J > 1) in Fig. 0(b). 

The time-dependence of the logarithmic negativity, as well as of the probability p, 
is clearly different in the two cases. It is illustrative to consider the case of only four 
sites, (with the two particles initially located at sites 2 and 3). Then the initial state is 
|10)sdq, which is a separable state of the two spatially delocalized qubits. If interactions 
are absent, then already at early times populations of states lying in the delocalized qubit 
space can originate from two possible paths (starting at sites 2 or 3, respectively). It 
is this interference which leads to a fast generation of entanglement. If interactions are 
strong, then one of this paths is effectively suppressed at early times as each particle is 
confined to "its" qubit (this leads to the larger probability p for a successful projection). 
Entanglement in this case is (initially) only generated through the collisional phase shift 
(this effect can be used to implement a phase gate for spatially delocalized qubits [27]). 
which however is small. As a consequence, entanglement is smaller if interactions are 
large. 

In the presence of a ground state background with an average number n of particles 
per site, generating (entangled) spatially delocalized quantum bits from the evolution 
of two extra particles only can be done effectively in the Mott case. Here the definition 
of the basis can be modified as |0)sdq = \ni = n+l,7ij+i = n), |1)sdq = \ni = n,n i+ i = 
n + l). In the superfluid case, the on-site particle number fluctuations in the ground 
state require a projection onto the subspace of fixed number of particles per spatially 
delocalized qubit, which strongly reduces the efficiency of the scheme. 
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4. Conclusions 

Summarizing, we have analyzed the generation of entanglement between the outer 
extremes of a ID Bose-Hubbard chain by loading an extra particle on top of the ground 
state. We have investigated effects arising from direct competition between tunneling 
and interactions in the entanglement behavior. In some limiting cases, the bosonic 
propagation can be adequately described as continuous time quantum walks. As part 
of our analysis, we have discussed two conceptually different "computational bases" 
to quantify entanglement between particles which have no internal degrees of freedom, 
namely the occupation number basis and a spatially delocalized qubit basis. 
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